Geometric Phase and Non-Adiabatic Effects in an Electronic Harmonic Oscillator 
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Steering a quantum harmonic oscillator state along cyclic trajectories leads to a path-dependent 
geometric phase. Here we describe an experiment observing this geometric phase in an electronic 
harmonic oscillator. We use a superconducting qubit as a non-linear probe of the phase, otherwise 
unobservable due to the linearity of the oscillator. Our results demonstrate that the geometric phase 
is, for a variety of cyclic trajectories, proportional to the area enclosed in the quadrature plane. At 
the transition to the non-adiabatic regime, we study corrections to the phase and dephasing of the 
qubit caused by qubit-resonator entanglement. The demonstrated controllability makes our system 
a versatile tool to study adiabatic and non-adiabatic geometric phases in open quantum systems 
and to investigate the potential of geometric gates for quantum information processing. 



It is a well known result in quantum mechanics that 
an adiabatically changing Hamiltonian causes a system, 
prepared initially in a stationary state, to follow the in- 
stantaneous energy eigenstates. This adiabatic theorem 
[H determines the evolution of the state up to a complex 
phase factor. Starting in an energy eigenstate, a cyclic 
adiabatic change of the Hamiltonian drives the system 
along a closed path in the jspace of physically distinct 
states. As noted by Berry y, even when corrected for 
the familiar dynamical phase given by the energy of the 
system, the initial and final state after the cyclic evo- 
lution can differ by a geometric phase determined only 
by the path. Geometric effects have been experimen- 
tally observed in a variety of two-level systems such as 
single neutrons [1] , nuclear spins [1, [1| and superconduct- 
ing qubits 0, 0] ■ The geometric phase of another simple 
quantum system - a harmonic oscillator - has been used 
for entanglinggates using harmonic motional modes of 
trapped ions [8] but has not been a subject of detailed 
experimental investigations in other systems. 

The independence of the geometric phase on dynami- 
cal quantities makes it stable under certain types of fluc- 
tuations in the system parameters f^-[ll|, offering inter- 
esting possibilities for potential noise-resilient quantum 
gates. Although the properties of geometric phases in 
the presence of noise in more general settings are still 
under debate 



12l . Il3j . geometric gates have been 



pro- 



posed in various physical implementations 
noise-induced geometric dephasing has been studied ex- 



perimentally [a, |23|, [2; 



Here we describe an experiment measuring the geomet- 
ric phase of an adiabatically manipulated harmonic oscil- 
lator in an electronic superconducting circuit 25[ . In con- 
trast to anharmonic systems such as qubits, the linearity 
of the harmonic oscillator implies that the accumulated 
geometric phase is the same for all energy eigenstates and 
therefore cannot be measured simply by observing the 
phase difference between two adiabatically transported 
eigenstates in a superposition 0. Instead, we utilize a 
dispersive coupling between a qubit and the oscillator 
which introduces a shift of the resonator frequency de- 



pending on the qubit state. The difference between the 
geometric phases accumulated for the two oscillator fre- 
quencies is then measured usin g th e qubit as an inter- 
ferometer, as also proposed in [26|. The good control- 
lability of our system allows us to investigate the phase 
for a wide range of parameters. In particular, we iden- 
tify corrections to the phase and dephasing of the qubit 
caused by qubit-resonator entanglement at the transition 
to non-adiabatic evolution. In this regime, fast geomet- 
rically protected gates based on the Aharonov-Anandan 
phase [271 . which are not restricted by the adiabaticity 
condition, may be realizable [28 1. 

In our setup, the harmonic oscillator is implemented 
as one of the electromagnetic modes of a transmission 
line resonator at a frequency of a;r/27r « 7.0 GHz (with 
the qubit in the ground state). It is dispersively coupled 
to a superconducting qubit of the transmon type (29| 
with an energy separation between the two lowest lying 
energy levels of hujq « /i x 8.3 GHz, an anharmonicity 
Q;/27r —0.4 GHz and a resonant coupling strength to 
the oscillator of .g/27r « 56 MHz. The sample is oper- 
ated in a dilution refrigerator at a base temperature of 
approximately 10 mK. 

We manipulate the state of the resonator using a mi- 
crowave drive field. The Hamiltonian in the reference 
frame rotating at the drive frequency w is 
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1J-|22[ and H = h{LJr - oj)a'^ a + -hei{t){a^ + a) + -hieQit){a'^ - a), 



where ej{t) and Sqit) are the in-phase and quadrature 
components of the drive, which we can control individ- 
ually. These two components are directly related to the 
amplitude e and phase tp of the drive by £/ -I- ieg = ee^^ . 

If the variation of ei[t) and £g(t) is slow compared 
with the detuning 5 — w,. — w, the resonator field adi- 
abatically follows the ground state of the Hamiltonian. 
This situation is analogous to the oscillations of a me- 
chanical harmonic oscillator, such as a mass on a spring, 
following the changes in amplitude and phase of a driving 
force [Fig.jTJa)]. The ground state is a coherent state |a), 
that is, an eigenstate of the annihilation operator a with 
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eigenvalue a = —{ej + ieQ)/25. Classically, the absolute 
value and phase of a correspond to the amplitude and 
phase of the resonator field, respectively. Its square |ap 
is the mean number of photons n in the resonator. 

We can manipulate the coherent state of the resonator 
at will by changing the drive components £/ and sq. 
These are calibrated in terms of the resonator pho ton 
number using a qubit ac Stark shift measurement [30| • If 
we make its quadrature a trace a closed path in the com- 
plex plane, returning back to the vacuum state a — 0, 
the system acquires a dynamical phase proportional to 
the time integral of + Eq and a geomet ric p hase —2 A, 
where A is the area enclosed by the path [31]. 

The field quadrature a depends on the detuning 6 be- 
tween the drive and the resonator frequency through the 
Lorentzian response function of the resonator character- 
ized by its center frequency ojr and its width k/2i: « 
500 kHz [Fig. 1(b)]. The dressed resonator frequency is 
determined by the state of the dispersively coupled qubit 
[ist . We denote it by for the qubit in its ground state 
\g). If the qubit is in the excited state |e), the resonator 
frequency is shifted to ujr + 2%. The size of the path 
traced by the coherent state, and hence also the accumu- 
lated phase 7'^-', therefore depends on the qubit state \s) 
[Fig. He)]. 

By applying a 7r/2 pulse to the qubit |32|, we initially 
prepare the system in the equal superposition of |0) (X) |g) 
and |0) €5 |e), where |0) is the resonator vacuum state. 
The adiabatic cycle takes the system into the state 
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(exp(i7(f))|0) ® \g) -t- exp(i7'"))|0) ® \e)). 



Note that due to the cyclicity and adiabaticity of the 
process, the resonator returns back to the vacuum state 
regardless of the state of the qubit and the two systems 
are again disentangled. We then finish the manipulation 
(Ramsey sequence; Fig. [2][a)) with a second 7r/2 pulse 
applied to the qubit, either in phase with the first one or 
shifted by 7r/2, and use dispersive readout [33| to measure 
its excited state population. In this way we obtain the x 
and y projections of the qubit state Bloch vector (ct), al- 
lowing us to calculate the phase difference 7 = 7^'^) — 7^^' 
as the angle of rotation of the Bloch vector about the 
z axis. This phase difference again contains a dynami- 
cal part, related to the ac Stark shift, and a geometric 
contribution 7^ = — 2AA, where is the area enclosed 
between the coherent state paths for the qubit in the 
ground and in the excited state [Fig. HJc)]. Drive am- 
plitudes typically used in our experiment result in up to 
n « 20 resonator photons, corresponding to 7g on the 
order of 27r. 

We control the resonator state to trace out different 
path shapes [Fig.lJJb)] in clockwise and counterclockwise 
direction and measure the accumulated phase. For each 
of these paths we also measure the phase when leaving 
the amplitude modulation unchanged, but keeping the 
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FIG. 1. (a) Adiabatically driven mechanical oscillator 
formed by a mass on a spring. Displacement is proportional 
to the quickly oscillating force with slowly varying amplitude 
and phase, (b) The Lorentzian resonator response function 
centered at frequencies ujr and LJr + 2x for the qubit in \g) 
and je), respectively, (c) Area of the coherent state path in 
the complex quadrature plane depending on the qubit state. 
The measured geometric phase difference is proportional to 
the area /S.A between the paths. Dashed circles with an area 
Avne = n/2 represent the size (rms) of vacuum fluctuations. 



phase of the drive constant [Fig. Hfc)]. The resulting 
straight paths yield the same dynamical phase as each of 
the original paths but no geometric phase since the area 
enclosed by them vanishes. 

In this way, the phases are measured for different 
lengths T of the drive pulse. In the adiabatic limit 
(T > 100 ns), the dynamical phase, measured for the 
straight path, scales linearly with T [Fig. [31[a)], as ex- 
pected. The geometric phase, evaluated as the differ- 
ence between the phase for the 'area-enclosing' path and 
the dynamical phase, approaches a constant value whose 
sign depends on the path orientation [Fig. [3Jb)]. As the 
cycle is traversed faster, and T becomes comparable to 
27r/|5| = 25 ns, non-adiabatic effects become apparent. 
In the weakly non-adiabatic regime (T > 25 ns), these 
are well described by an expansion in powers of 1/T [33 |. 
Notably, the observed non-adiabatic corrections remain 
small even if the adiabaticity condition T ^ 27r/|(5| is 
clearly violated. 

The scaling of 7^ with the area AA (determined by 
the amplitude of the drive pulse) is observed to be lin- 



3 



-a 

"D- 

£ 



Qubit drive 



-I — ih 



Resonator drive 



Non-adiabatic 



Adiabatic 



Readout 



Sq 




8 


|A. 






S| 


t 




FIG. 2. (a) Illustration of the qubit and resonator drive 
pulses. The maximum amplitude of the resonator drive is de- 
noted by eo. (b) The four different shapes of the path used in 
the experiment and the corresponding geometric phases rela- 
tive to the one of the circular path, (c) Path of the resonator 
drive components in the IQ plane (left) and the corresponding 
time dependence (right) of the components ei and eq (solid 
line) and the drive amplitude (ef+eg)^'^^ (dashed line) for the 
circular drive path (top row) and the straight path (bottom 
row) yielding the same dynamical phase. 



ear [Fig. |31Ic)], in agreement with theory. We also ver- 
ify this scaling by comparing the adiabatic geometric 
phase, obtained by extrapolating the measured geomet- 
ric phases to T — > oo, for different path shapes indicated 
in Fig.[2][b). The ratios of the extracted phases to that of 
the circular path are 0.493±0.016 (1/2) for the semicircle, 
0.647±0.016 (2/7r w 0.637) for the square and O.OOiO.07 
(0) for the figure-eight shape, in excellent agreement with 
the theoretical values stated in parentheses. 

We have also measured the dependence of the geomet- 
ric phase on the drive detuning 5, observing its increase 
with decreasing 5. This trend is explained by the larger 
field amplitudes and, thus, an increased geometric phase 
when driving the resonator closer to its resonance fre- 
quency. Changing the sign of the detuning 5 reverses the 
sign of the resulting phase since the relative size of the 
amplitudes of the resonator field for the two qubit states 
is reversed. The measured dependence agrees well with 
these theoretical predictions [Fig. EJd)] . 

Apart from the rotation of the qubit state Bloch vec- 
tor about the z axis, representing the accumulated phase 
difference, we also observe a decrease in the length R 
of its xy projection, i.e. dephasing of the qubit. The 
measured value of R shows a strong dependence both 
on the drive pulse amplitude [Fig. SJa)] and duration 
[Fig. m^b)]. The dephasing effect can be explained as 
a result of non-adiabaticity which leads to entanglement 
between the oscillator and the qubit. For a non-adiabatic 
drive cycle, the resonator coherent state does not fol- 
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FIG. 3. (a) Measured total phase difference 7 for a res- 
onator pulse of fixed maximum amplitude eo/Svr « 370 MHz, 
detuning 5/2tt = 40 MHz and varying duration T, tracing 
a counterclockwise circular (green circles), clockwise circu- 
lar (orange squares) and straight trajectory (blue diamonds) 
with the same time-dependence of the drive amplitude. Solid 
lines show theory, (b) Dependence of the geometric phase 
on the drive pulse duration for the counterclockwise and 
clockwise trajectory. Solid lines show the adiabatic limit, 
dashed lines are fits with corrections proportional to l/T and 
(c) Geometric phase for a counterclockwise and clock- 
wise circular pulse with fixed duration T — 300 ns, detuning 
i5/27r = 40 MHz and varying maximum amplitude, plotted 
as a function of the area AA enclosed between the coherent 
state trajectories, shown together with fitted linear functions 
(solid lines). Dashed lines represent the theoretical depen- 
dence 7g = —2AA. (d) Measured geometric phase in the adi- 
abatic limit as a function of the detuning 5, compared with 
the theoretical dependence (solid line). 



low the changes of the drive parameters [Fig. HJc)] and 
its trajectory exhibits periodic excursions from the adia- 
batic path, i.e., ringing. The system does not necessarily 
return to a product state of the qubit and the resonator 
and in general, the final quadratures of the resonator field 
for the two qubit states differ. Therefore, the two sub- 
systems remain entangled 



and the reduced density 
matrix of the qubit corresponds to a mixed state with 
a Bloch vector length decreased by a factor given by the 
overlap 
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of the two final resonator states. Alternatively, this 
dephasing effect can be understood as measurement- 
induced dephasing - a result of the resonator field ex- 
tracting information about the qubit state (36j . 

The Gaussian dependence of R on the drive amplitude 
[Fig. mja)] follows directly from Eq. (1) and the propor- 
tionality between a and e. The fall-off of R is faster 
for shorter, less adiabatic pulses as they result in larger 
separation between the final field quadratures and hence 
stronger dephasing. Interestingly, the magnitude of the 
dephasing effect shows oscillations in the evolution time 
T and also depends on the orientation of the drive cy- 
cle, as seen in Fig. |4l[b). The asymmetric behavior has 
a simple explanation in terms of the Fourier transform of 
the drive signal. For the circular path, the two dominant 
Fourier components have frequencies w and uj ± 27r/T, 
where the signs + and — correspond to the clockwise and 
counterclockwise orientation of the path, respectively. In 
one of these two cases, the resonator is driven closer to 
its resonance frequency, resulting in stronger dephasing. 
For fast clockwise drive cycles (T ^ 50 ns) the qubit 
is almost completely dephased and the projected Bloch 
vector length R approaches zero. The oscillatory be- 
haviour of R with a period of 25 ns corresponding to 
S/27T — 40 MHz, is due to the ringing of the resonator 
field. For cycle periods corresponding to the maxima 
of the Bloch vector length R, at which the ringing fre- 
quencies are commensurate, the resonator state path is 
cyclic and dephasing is minimized. These particular peri- 
ods could be used to realize fast non-adiabatic geometric 
gates. 

The high level of control over the coherent state of the 
resonator field has allowed us to measure its geometric 
phase for a wide range of path shapes in the adiabatic 
and the non-adiabatic regime. The characteristic features 
of the geometric phase that we have observed - indepen- 
dence on dynamical quantities and scaling with the area 
enclosed by the trajectory in the parameter space - are in 
good agreement with theory. Non-adiabatic effects intro- 
duce corrections to the geometric phase and dephasing of 
the reduced qubit state due to residual qubit-resonator 
entanglement. In the future, our system can serve as 
a valuable tool for investigations of non-adiabatic geo- 
metric phases to shed light on their robustness and po- 
tential use in fast geometric quantum gates. 
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